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Abstract. It is well-known that orthogonal polynomials on the real line satisfy 
a three-term recurrence relation and conversely every system of polynomials satis- 
fying a three-term recurrence relation is orthogonal with respect to some positive 
Borel measure on the real line. In this paper we extend this result and show that 
every system of polynomials satisfying some (2A'' + l)-term recurrence relation can 
be expressed in terms of orthonormal matrix polynomials for which the coefficients 
are N x N matrices. We apply this result to polynomials orthogonal with respect 
to a discrete Sobolev inner product and other inner products in the linear space of 
polynomials. As an application we give a short proof of Krein's characterization 
of orthogonal polynomials with a spectrum having a finite number of accumulation 
points. 



1. Introduction 

A sequence of orthonormal polynomials Pnix) (n = 0, 1, 2, . . . ) on the real line, 
orthonormal with some probability measure /x, always satisfies a three-term recur- 
rence relation 

(1.1) XPn(x) = ttn+lPn+lix) + b^Pnix) + a„p„_i(x), 

with initial conditions p-i{x) = and po{x) = 1. The recurrence coefficients are 
given by 

The converse is also true: a system of polynomials satisfying a three-term recurrence 
relation (1.1) with ttn+i > and 6„ G M (n = 0, 1, 2, . . . ) is always a system of 
orthonormal polynomials with respect to some probability measure /i on the real 
line. This converse result was given by Favard in 1935 [10] but was known earlier 
and appears already in the books by Stone [21, Theorem 10.27 on pp. 545-546], 
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Perron [18, §36, Satz 4.6/4.7] and Wintner [22, §32 and §87] who attributes the 
case of finite support to E. Heine [12, §108]. In an elementary form (lacking the 
Riesz representation theorem) the result is in Stieltjes' 1894 work [20, §11] and 
apparently the result was already known by Chebyshev [5] when the support of the 
measure /x is finite. 

The importance of Favard's theorem is that orthogonal polynomials and poly- 
nomials satisfying a three-term recurrence relation are the same thing. Properties 
regarding zeros (real and simple zeros, interlacing of zeros) and positivity of con- 
nection coefficients can thus be studied from two points of view: on one hand using 
the orthogonality, on the other hand using the recurrence relation. 

Favard's theorem can also be proved for orthogonal matrix polynomials. Or- 
thogonal matrix polynomials on the real line have been considered in detail by 
M.G. Krein [13]. See also the book by Berezanskii [4] and more recent papers 
by Aptekarev and Nikishin [1], Geronimo [11], and Sinap and Van Assche [19]. 
Consider matrix polynomials satisfying the three-term recurrence relation 

(1.2) xPkix) = Dk+iPk+i{x) + EkPkix) + DlPk-i{x), 

with Pq{x) = I and P-i{x) — 0, where Pk{x) are matrix polynomials with coeffi- 
cients in C^^^ and the recurrence coefficients Dk+i-,Ek are also NxN matrices for 
which = Ek and dct 7^ 0. By using spectral theory Aptekarev and Nikishin 
[1] show that the polynomials Pn{x) (n = 0, 1, 2, . . . ) are orthonormal with respect 
to some Hermitian matrix of measures M = {l-i'k,i)k 1=1 which is positive definite: 

J P^{x) dM{x)P:;^{x)^5^,nI- 

Recently one of us studied polynomials Pn{x) satisfying a {2N + l)-term recur- 
rence relation 

N 

(1.3) h{x)pn{x) = Cn,QPn{x) + ( ^n,kPn—k (^n+k,kPn+k 

fc=l 

where /i is a polynomial of degree and Cn^k (n = 0, 1, 2, . . . ) are real sequences for 
k — 0,1, . . . ,N with Cn^N 7^ 0. From [7] we get, after straightforward reformulation, 
that a sequence of polynomials satisfies a (2iV -|- l)-term recurrence relation if and 
only if the following orthogonality condition holds: there exists a N x N matrix of 
measures = {f^k,l)kr=o such that the bilinear form 

N-l 

(1-4) Bi^{p,q)=^ I Rh,N,k{p)Rh,N,l{q)dlik,h 

k,l=0'^ 

where 

m N—1 m 

(1.5) Rh,N,kiP)i^) = J2^'''^^'' p{x) =^^ak,nx''h'^ix) 

n=0 fc=0 n—O 

is an inner product on the linear space of polynomials P and {pn)n is the sequence 
of orthonormal polynomials with respect to B^. In [8] this is improved showing 
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In Section 2 we will show that polynomials satisfying a higher order recurrence 
relation of the form (1.3) (but with complex coefficients) are closely related to 
matrix polynomials satisfying a thrcc-tcrm recurrence relation and that Favard's 
theorem for matrix polynomials and Favard's theorem for polynomials satisfying a 
higher order recurrence relation are the same. 

In Section 3 we show how discrete Sobolev orthogonal polynomials and some 
new orthogonal polynomials with respect to an inner product of the form 

(/, ^) = J fi^)9ix) dn{x) + o,kf{ck)^ afe5'(cfe) j , Ofe, Cfe e IR 

are related to orthogonal matrix polynomials and we explicitly give the orthogo- 
nality matrix of measures M in terms of the parameters in the inner product. The 
advantage of working with the matrix polynomials is that the orthogonality con- 
ditions no longer require the evaluation of a function and its derivatives in various 
points. Finally, in Section 4 we show how Krein's theorem regarding orthogonal 
polynomials with a spectrum with finitely many accumulation points follows easily 
by this correspondence between higher order recurrence relations and orthogonal 
matrix polynomials on the real line. 

The connection between matrix polynomials and scalar polynomials obtained 
in the present paper shows some analogy with the connection between orthogonal 
polynomials on a lemniscate and orthogonal polynomials on the unit circle, as 
given by Marcellan and Rodriguez [15] and the connection between orthogonal 
polynomials on an algebraic harmonic curve and orthogonal matrix polynomials on 
the real line, as given by Marcellan and Sansigre [17]. In fact the particular choice 
of basis in the linear space of polynomials that we use in Sections 2 and 3 is the 
same as the basis used in [15] and [17], but the polynomial h{x) in the present paper 
is connected with the recurrence relation, whereas in [15] this polynomial describes 
the lemniscate and in [17] it describes the algebraic harmonic curve. 

2. Recurrence relation and matrix polynomials 

In order to establish the main theorem we need to consider the operators Rh,N,k 
defined in the introduction of this paper (see (1-5)). For the sake of simplicity we 
start by considering the case h{x) = . Then, we denote the operators by RN,m 
and it is not hard to see that they are defined by 

i.e., the operator RN,m takes from p just those powers with remainder m modulo 
N and then, removes x'^ and changes x^ to x. Thus, we have 

p{x) = i?Ar,o(p)(x^) + xi?iv,i(p)(x^) + • • • + x''-^Rn,n-i{p){x''). 

For example, when N = 3 and p{x) = 6x^ + 5a;^ + Ax^ + Sx'^ + 2a; + 1 we have 
R3,oip)i^) = 4:X + 1, Rs,i{p){x) = 5x + 2, Rz^2{p){x) = 6x + 3. For iV = 2 this 
corresponds to taking the odd (i?2,i) and even (-^2,0) parts of the polynomial p. 
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Theorem. Suppose Pn{x) (n — 0,1,2, ■ ■ ■ ) is a sequence of polynomials satisfying 
the following (27V + l)-term recurrence relation 

N 

(2.1) X^Pn{x) = Cn,0Pn{x) + ^ (c^Pn-fe(a^) + Cn+k,kPn+k{x)) , 

k=l 



where Cn,o {n — 0,1,2, . . .) is a real sequence and Cn,k (n = 0, 1, 2, . . . ) are complex 
sequences for k = 1, . . . , N with Cn^N 7^ for every n and with the initial conditions 
Pk{x) = for k < and pk given polynomials of degree k, for k = 0, ■ ■ ■ , N — 1. 
We define the sequence of matrix polynomials {Pn)n by 



Pn{x) = 



/ RN,o{PnN){x) ••• RN,N-l{PnN){x) \ 

RN,o{PnN+l){x) ■■■ RN,N-l{PnN+l){x) 

^RN,oiPnN+N-l){x) ■■■ RN,N-l{PnN+N-l){x) / 



Then this sequence of matrix polynomials is orthonormal on the real line with respect 
to a positive definite matrix of measures and satisfies a matrix three-term recurrence 
relation. Conversely, suppose Pn = {Pn,m,j)^~j=o ^ sequence of orthonormal 
matrix polynomials or equivalently satisfying a matrix three-term recurrence relation 
(without loss of generality we can assume the leading coefficient of Pn to be a lower 
triangular matrix), then the scalar polynomials defined by 

N-l 

(2.2) PnN+m{x) = J2 ^'Pn,mA^'')^ (n G N, < m < AT - 1), 
satisfy a {2N + l)-term recurrence relation of the form (2.1). 

Proof. The equivalence between {pn)n satisfying a (2A^+l)-term recurrence relation 
and {Pn)n being a sequence of matrix orthonormal polynomials is a consequence 
of the definition of the Pn from the pk (or conversely, the pk from the Pn) and the 
orthogonality condition (1.4). 

Let us show that the matrix polynomials (Pn)n satisfy a matrix three-term re- 
currence relation, giving explicitly the matrix coefficients which appear in this re- 
currence formula. To do that, we consider the A'"-Jacobi matrix J associated to 
{Pn)n, which is the (2A^ + l)-banded infinite Hermitian matrix defined by putting 
the sequences (c„ ;)„ which appear in the recurrence relation on the diagonals of 
the matrix J, i.e., we define the matrix J = {jn,m)n,men by 

if < n - m < A^, 



-'n, n— m 



jn,m = { Cm,\n-m\ iiO<m-n<N, 

if |n - m| > N. 

Now, we split up this A^-Jacobi matrix in blocks of dimension N x N, and then we 
get the N X N matrices En and Dn defined by 

. • / CnN+i,\i-l\ ifi>/, 

[-t^n id — 3nN+i,nN+l — \ .n ■ ^ , n >U, 
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and 



{Dn)i,l = 



ifi</, 

j{n-l)N+i,nN+l if ^ > /, 

ifi</, 

CnN+l,N+l-i Hi > I, 



n > 1, 



I.e., 



with 



J 



/ CnN,N 
CnN,N-l 
CnN,N-2 



El 



D2 

E2 



\ 





CnN+l,N 
CnN+l,N-l 






CnN+2,N 



CnN,l CnN+1,2 CnN+2,3 









CnN+N-1, 



1,N ' 



and 



En — 



CnAr+2,2 



CnAf+1,0 
Cn Ar+2,1 



CnAf+2,2 
Cn Ar+2,1 
CniV+2,0 



CnAr+Af-l,Af-l \ 
CnAr+iV-l,Af-2 
Cn Ar+Ar-l,Af-3 



CnTV+iV-l.JV-l CnAf+Ar-l,Ar-2 CnAr+Ar-l,Af-3 ••• CnAr+Ar-1,0 ' 

Then E'^ is an Hermitian matrix, and the conditions c^^jv 7^ shows that is a 
lower triangular matrix with det ^ 0. 
If we compute the following expression 

(2.3) L>n+iPn+i(x) + EnPn{x) + DlPn-x{x) 

we find that the entry (A;, m) (0 < A;, m < iV — 1) of this matrix is equal to 

k 

1 ) jv_|_ j ^ jv_|, j _ fc -Rjv, nx (P(n+ 1 ) iV+ j ) (^) 
j=0 j=k 



3 

N-1 



+ XI CnN+k,N+k-jRN,m{P{ n—l)N+j 

that is, 

(2.4) CnN+k,oRN,m{PnN+k){x) 
N 

+ X {CnN+k,lRN,m (PnN+k-l ) (a^) + CnJV+A;+/,«-RiV,m(PniV+fc+/) (^)) • 
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But from the (2iV+l)-term recurrence relation which the polynomials {pn)n satisfy, 
it follows that the sequence of polynomials {RN,i{Pn))n = 0, • • • ,N — 1) satisfies 
the following recurrence formula 

N 

xRN,i{Pn){x) = C„fiRN,i{Pn){x) + [c^iRn ,i{Pn-l){x) + Cn+l,lRN ,i{Pn+l){x)) 

1=1 

and so (2.4) is equal to xRN^rn{PnN+k){x). Hence (2.3) is equal to xPn{x) and we 
have proved that the sequence of matrix polynomials satisfies the matrix three-term 
recurrence relation 

Note that the matrix polynomials obtained in this way always have a leading coef- 
ficient which is a lower triangular matrix and also the matrices Dn which appear 
in the recurrence formula are lower triangular. 

To prove the converse, suppose {Pn)n is a sequence of matrix polynomials satis- 
fying the recurrence formula 

xPn{x) = Dn+lPn+l{x) + EnPn{x) + D*^Pn-l{x), 

where det ^ 0. First, we prove that we can assume both the leading coefficients 
of Pn and the matrices to be lower triangular matrices, which is necessary since 
otherwise the scalar polynomials Pn defined by (2.2) could have degree different 
from n. Orthonormal matrix polynomials for a given orthogonality measure M 
are only determined up to a unitary factor in the sense that UnPn{x) is also 
orthonormal with respect to the measure M whenever UnU* — I. If Qn{x) is a 
sequence of orthonormal matrix polynomials satisfying 

xQn{x) = An+iQn+i{x) + BnQn{x) + A*Q„_i(a;), 

then the polynomials Pn{x) — UnQn{x), with UnU* — I, satisfy 

xPnix) = UnAn+lU*+^Pn+iix) + UnBJJ^Pnix) + t/„<C/*_iP„_i (x) . 

If the non-singular matrices are not lower triangular, one can always find unitary 
matrices U^ such that 

Dn — Un—lArJJ^ 

are lower triangular matrices. These unitary matrices are given recursively as fol- 
lows. First we choose a unitary matrix Uq such that UqQq = Pq is lower triangular. 
Next we observe that the matrix UqAi can always be factorized as UqAi = DiUi, 
where Ui is a unitary matrix and Di is lower triangular (QR-factorization of Francis 
and Kublanovskaja). With this choice of Ui we thus have Di = UqAiUi- In gen- 
eral, if Uq, Ui, . . . , Un-i have been obtained, then we apply the QR-factorization to 
Un-iAn to find a unitary matrix U^ with Un-iA^ = D^Un, with lower triangu- 
lar. Therefore any system of orthonormal matrix polynomials can be transformed 
to a system of orthonormal matrix polynomials with leading coefficients which are 
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Now, since we can assume the leading coeflBcient of P„ to be lower triangular 
with non- vanishing determinant, the sequence of scalar polynomials {pn)n defined 
by (2.2) are such that dgr{pn) = n. And proceeding as before, it is no hard to prove 
that they satisfy a (2A^ + l)-term recurrence relation like (2.1). The coefficients in 
this recurrence formula can be obtained from the matrices D^, E^. □ 

For the general case of a (2A'" + l)-term recurrence relation defined by a polyno- 
mial h of degree N 

N 

h{x)pn{x) = CnfiPn{x) + ^ ( Cn,kPn—k (^n+k,kPn+k 

fc=l 

the theorem works again if we change the operators RN,m to Rh,N,m (m = 0, • • • , N— 
1) (see (1.5)) in the definition of the matrix polynomials from the scalar polyno- 
mials pn- Thus, in the definition of the operators Rh,N,m, instead of using the basis 
of monomials {1 ^, . . . } to span the linear space of polynomials, we will use 

the basis 

(2.5) {l,a;, ... ,x^~^,h{x),xh{x),... ,x^~^h{x),h^{x),xh'^{x), . . .} 

= {x'^h'^ix) : n = 0, 1, . . . , - 1, m = 0, 1, 2, . . . }. 

A polynomial p of degree nN + m (0 < m < N) can then be expanded in this basis 
as 

n N-1 
i=0 j=0 

Now, the operator Rh,N,j takes from p just those terms of the form aijx^h^{x) and 
then, removes the common factor x^ and changes h{x) to x. 
Conversely, since now we have 

p{x) = Rh,N,o{p){h{x)) + xRh,N,i{Kx)) H h x^~^Rh,N,N-i{p){Kx)) 

we change (2.2) in the theorem to 

N-l 

PnN+m{x) = Yl X^Pn,m,j{h{x)), (n G N, < m < AT - 1) 

in the definition of the scalar polynomials pn from the matrix polynomials Pn- 

3. Examples 

3.1 Discrete Sobolev orthogonal polynomials. 

An important class of polynomials satisfying a higher order recurrence relation 
is obtained by taking polynomials orthogonal with respect to an inner product of 
(discrete) Sobolev type 

N 

{f,9)= / /(a;)^(^)M^) + 5]Afc/(^'=)(cfe)g(^'=)(cfc), 
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where are positive real numbers and are real numbers (which are allowed to 
coincide) . For a discrete Sobolev inner product 

{f.9) = J f{x)g{x)dpi{x) + \f^-\c)g^'\c) 

it was shown by MarccUan and Ronveaux [16] that the corresponding monic or- 
thogonal polynomials qn{x) satisfy a (2r + 3)-term recurrence relation 



n+r+l 

{x-cY+'^qn{x) = InjQjix). 

j=n—r—l 



The orthonormal polynomials then satisfy a recurrence relation of the form (1.3) 
with h{x) = {x — cY'^^ and thus these polynomials are also orthogonal with respect 
to a positive definite matrix of measures. Evans et al. [9] consider weighted Sobolev 
inner products 

(/'^) = E / f'\^)9^'\^)dfik{x), 

k=0'' 

and show that the existence of a polynomial h for which 



(3.1) {hp,q) ^ {p,hq) 

for all polynomials p, q implies that the measures fik < k < N) are necessarily 
discrete with support at the zeros of h and that the orthonormal polynomials satisfy 
a (2A^+ l)-term recurrence relation of the form (1.3), where N is the degree of the 
minimal polynomial h for which (3.1) holds. It follows that polynomials which are 
orthonormal with respect to a discrete Sobolev inner product correspond to matrix 
polynomials orthonormal with respect to a positive definite matrix of measures. 
Consider the discrete Sobolev inner product 



/M Mi 
p{x)q{x) dfi{x) + KjP^'Hci)q^'\ci), 
i=i j=i 

where p, q are polynomials, Xi^j > and N — M + Mj. Here derivatives are 

taken at M points Cj e IR and at the point Cj the highest derivative is of order Mj. 
Introduce the polynomial 



M 



H^) Yli^ - Ci) 



Mi + 1 



then h is of degree and has its zeros at the points q where the derivatives of the in- 
ner product are evaluated. Instead of using the basis of monomials {1, a;, a;^, x^, . . . } 
to span the linear space of polynomials, we will use the basis 

(3.3) {l,a;, ... ,x^~^ ,h[x)^xh{x), . . . ,x^~^h{x),h^{x),xh^{x), . . .} 
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A polynomial p of degree Nk + 1 {0 < I < N) can then be expanded in this basis as 

N-l k 



n—0 m—0 

where a„^fc = whenever n > I. Taking the terms in x'^ together and putting 

k 



Rh,N,n{p){x) = an, 



m=0 



then for < n < iV each polynomial Rh,N,n{p) has degree at most k (for n > I the 
degree is less than k) and 



N-l 



p{x) = J2 x''Rh,N,n{P)iKx))- 



n=0 



The polynomial p is thus equivalent (modulo h) with the vector polynomial given 
by {Rh,N,o{p): Rh,N,iip), ■ ■ ■ , Rh,N,N-iip)) and we will write 

P = {Rh,N,o{P)^ Rh,N,l{p), ■ • ■ , Rh,N,N-l{p))- 

Observe again that in case h{x) = x^ this amounts to the decomposition of a 
polynomial p into its odd and even parts. 

For the jth derivative of p we then find, using Leibnitz' rule for the derivative 
of a product, 

^ /A\ M-k 

pU)(x) = ^ ^ P jn(n - 1) ■ ■ ■ (n - A; + l)x-'^ ^—^R^^^^^{p){hix)). 

n=j fc=0 ^ ^ 

The derivative of the composite function Rh,N,n{p){h{x)) can be evaluated using 
the formula of Faa di Bruno, giving 

(I'm 

Rh,N,n{p){h{^)) 



dx"^ 

m 



where are non-negative integers satisfying ai + a2 + ■ ■ ■ + am = i and ai + 2a2 + 
• • • + mam = m. If m < Mi and if we evaluate this expression at Cj then due to the 
fact that h'{ci) = • • • = h^^\ci) = we see that 



Rh,N,n{p){h{x)) 



dx 



0, 1 < m < Mj 



Therefore the only contribution in the expression for p^^\ci) when 1 < j < Mj is 
when k = j, giving 



N-l 

n 



p^'\ci) = E rTr^^r''^/.,iv,n(p)(o), 1 < i < M,. 
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For p = (i?/,,Ar,ob), • ■ • ,Rh,N,N-iip)) and q = (i?ft,7v,o(g), ■ • ■ , Rh,N,N-iiq)) the 
inner product (3.2) can thus be written as 



{p,q) = J {Rh,N,o{p)iKx)) ■■■ Rh,N,N-i{p)ih{x)))dM{x) 



/ Rh,N,o{(l)ihix)) \ 

Rh,N,i{(i)iK^)) 
\Rh,N,N-i{q){h{x)) / 



+ {Rh,N,o{PM ■■■ Rh,N,N-l{p){0))l^ 



where M is the N x N matrix of measures 



/ ^ft,iv,o(9)(0) \ 
Rh,N,i{qm ^ 

\Rh,N,N-iiqM^ 



(3.4) 



dM{x) = 



I d\x{x) X d^{x) 

X dn{x) x^ dn{x) 

x^ diJ,{x) x^ dji^x) 

\ x^~^ diJ,{x) x^ dii{x) 



X 



N-l 



dfx{x) \ 
x^ dfi{x) 
x^~^^ dji^x) 



X 



2N-2 



dfi{x) / 



and L is the matrix 



M Mi 



with L{i,j) the N x N matrix 
/ \ 







{k-jy. 



(0...0 



fc! k-j 



{k-jy.^i 



(N-iy N-i-j \ 
{N-i-jyH ) 



(AT-l)! N-l-j , 
\(iV-l-i)!S / 



with entries 



Lfc,n(«,j) 



k\n\ 



ji+k-2j 



{k-j)\{n-j)\ 



j <k,n< N -1. 



It is clear that L(i, j) is positive definite, hence L is also positive definite when all 
K,j > 0. If Pn{x) {n = 0, 1,2, . . .) are the orthonormal polynomials with respect 
to the Sobolev inner product (3.2), then we can write the polynomials pkN+i{x) 
< I < N using the basis functions (3.3): 



JV-l 



PkN+l{x)= x''Rh,N,n{PkN+l){h{x)), I = 0, 1, . . . , N - 1. 
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The matrix polynomials 



Pnix) 



/ Rh,N,o{PnN){x) ■■■ Rh,N,N-l{PnN){x) \ 

Rh,N,o{PnN+l){x) ■■■ Rh,N,N-l{PnN+l){x) 

\ Rh,N,o{PnN+N-l){x) ■■■ Rh,N,N-l{PnN+N-l){x) / 



are then orthonormal with respect to the matrix of measures M(/i~^) to which 
a mass point at is added, with weight given by the matrix L. The matrix of 
measures M(/i~^) is given by 



fW<<M(,r.M)^/f,M.))dM,.). 



where F : R — > 'R^ is a vector function such that F{h) e Li(M). In this way 
Sobolev orthogonal polynomials (with a discrete Sobolev part) can always be ex- 
pressed as orthogonal matrix polynomials, where the spectral matrix of measures 
has a mass point at the origin. Note that the matrix polynomial Pk{x) is of degree 
k with a leading coefficient which is a lower triangular matrix because the degree 
of Rh,N,n{PkN+i) is less than k whenever n > I. 

3.2 Perturbation of a measure on the real line by finitely many function 
values. 

As a second example we consider the inner product 

(3.5) {f,g) = / f{x)g{x)dfx{x)+ I ^ akf{ck) j akg{ck) J , 

\k^o J \k=0 J 

where /, g are real functions, is a positive measure on the real line and afc, Cfc G M 
{k = 0, . . . , N — 1). When = 1 then this is just a an inner product in L2 
with respect to the measure /i to which a mass Gq is added at the point cq. The 
case N = 2 with gq — —ai is basically an inner product involving differences, as 
introduced recently by Bavinck [2] [3] . 
If we consider the polynomial 

N-l 
fc=0 

then h vanishes at the points Cfc where the function values are evaluated in the 
inner product, and hence 

{f,hg) ^ J h{x)f{x)g{x)dfi{x) = {hf,g), 

which immediately implies that the orthonormal polynomials {pn)n for this inner 
product satisfy a (2^" + l)-term recurrence formula of the form (1.3). Expanding 
a polynomial p in the basis (3.3) gives 

JV-l 

Pi^) = x''Rh,N,n{p){h{x)). 
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Observe that 



N-l 



n=0 



hence for p = {Rh,N,o{p)^ ■ • ■ > Rh,N,N-i{p)) and q = {Rh,N,o{Q)^ • ■ • > Rh,N,N-i{q)) 
we have 



'AT-l 



Y akf{ck) Y akg{ck) 



, /c=0 



= (i?h,iv,o(p)(0) 



i??.,iv,iv-i(p)(0))L 



/ Rh,N,o{Q)iQ) \ 
\Rh,N,N-i{q){0)' 



where 



lZ^fc=0 ^fe 2^k=0 ^kCk ■■■ Z^fc=0 ^feCfe J- 



AT-l 



The inner product (3.5) in the hnear space of scalar polynomials thus corresponds 
to an inner product in L2 for the measure Nl{h~^), with M given by (3.4), to which 
a Dirac measure at is added with weight given by the positive definite matrix L. 
Moreover, the matrix polynomials 



Pnix) 



I Rh,N,o{PnN){x) 
Rh,N,o{PnN+l){x) 



Rh,N,N-l{PnN){x) \ 
Rh,N,N-l{PnN+l){x) 



^Rh,N,o{PnN+N-l){x) ■■■ Rh,N,N-l{PnN+N-l){x) J 

are orthonormal with respect to M(/?,~^) + IiSq. 

The special case when one adds to the i^2(A*) inner product a part dealing with 
the {N — l)st difference, i.e., 

iV-l 

J2 akfick) = A^-V(c), 



fc=0 



corresponds to the choice Ck = c + k5 and ak = j. )■ Observe that in this 

case 

/ \ 



: (0 ... {-5)^-^{N 



V(-<5)^-i(iV-l)!/ 

so that L contains zeros everywhere, except for the entry in the lower right corner. 
This case covers the inner products considered by Bavinck [2] [3]. The limiting 
case where 5 — > corresponds to a discrete Sobolev inner product with a {N — l)st 
derivative at the point c. The general discrete Sobolev inner product (3.2) can also 
be obtained as a limiting case of (3.5) by letting some of the Ck coincide and taking 
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4. Krein's theorem revisited 



Suppose that the spectrum cr(J) of a (tridiagonal) Jacobi matrix J is denu- 
merable. If the set of accumulation points of the spectrum a (J) is finite, then a 
complete characterization of the derived set (t(J)' is given by M.G. Krein [14] [6, 
Chapter IV, Section 6]: 

Krein's theorem. Suppose J is a hounded Jacobi matrix. Then every accumula- 
tion point of o-{J) is a zero of the polynomial h{x) of degree N if and only if the 
operator h{ J) is compact. 

The relationship between polynomials satisfying a higher order recurrence re- 
lation an orthogonal matrix polynomials enables us to give a short proof of this 
result. 

Proof of Krein's theorem. The infinite matrix h{J) is a banded matrix with band 
width 2N + 1, and can thus be written as 

N N 

where D, Ai, Bi, . . . , An, are diagonal matrices and V is the shift operator that 
acts as {y^l))n — '071+ 1- A diagonal matrix is compact if and only if its entries tend 
to zero, and compact operators on a Hilbert space form a two-sided ideal in the 
set of bounded operators on this Hilbert space. Hence the boundedness of the shift 
operator implies that h{ J) is compact if and only if the entries of h{J) tend to zero 
along the diagonals. If we write h{J) as a block Jacobi matrix 



h{J) = 



V 



El 



D2 

E2 



I 



where -Efc, are NxN matrices and are lower triangular, then the compactness 
of h{J) implies that Ef~ and D/. converge towards the zero matrix. A compact 
operator has a denumerable spectrum with zero as the only accumulation point. If 
/i is the spectral measure for J (with corresponding orthogonal polynomials Pn{x)), 
then M(/i~^) is the matrix of measures for which the matrix polynomials Pn{x) 
corresponding with h{J) are orthogonal, where 

dn{x) xdn{x) ■•■ x^~^ dfi{x) \ 



( 



dM(x) = 



xd[i[x) X d\i{x) 



X- 



diJ,{x) 



X' 



' d^{x) 



\x^ ^ dii{x) x^ d^x{x) 



X 



N 



X 



N+1 



d^{x) 
dfj,{x) 



X 



2N-2 



dfj.{x) / 



If a{h{J)) is the spectrum of h{J), then it follows that the spectrum of J is given 
by h~^{a{h{J)). We know that a{h{J)) has only one accumulation point at zero, 
which therefore corresponds to accumulation points of J at h~^{0), which are the 
zeros of the polynomial h{x). 

Conversely, if we know that cr(J) is denumerable with accumulation points at 
the zeros of h, then the spectrum of h{J) is given by h{a{J)) and is therefore 
also denumerable with accumulation points at h{h~^{0)) = 0, which makes h{J) a 



n 
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